Coherent control of optical four-wave mixing by two-color uj-3uj ultrashort laser pulses 
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A theoretical investigation on the quantum control of optical coherent four-wave mixing interac- 
tions in two-level systems driven by two intense synchronized femtosecond laser pulses of central 
angular frequencies lu and 3lu is reported. By numerically solving the full Maxwell-Bloch equa- 
tions beyond the slowly-varying envelope and rotating-wave approximations in the time domain, 
the nonlinear coupling to the optical field at frequency bui is found to depend critically on the 
initial relative phase cA of the two propagating pulses; the coupling is enhanced when the pulses 
interfere constructively in the center ((f) — 0), while it is nearly suppressed when they are out of 
phase ((f) — 7r). The tuning of the initial absolute phase of the different frequency components of 
synchronously propapagating ui-3ui femtosecond pulses can serve as a means to control coherent 
anti-Stokes Raman (CARS) processes. 



In recent years, encouraged by the developments in 
the engineering of intense ultrashort laser fields with a 
well defined absolute phase , studies on the phase con- 
trol of the interaction of two-color strong ultrashort laser 
pulses in nonlinear materials have received a great in- 
terest 0, S IH 0, 0- Phenomena arising from such 
ultrashort pulse interaction can be of extreme impor- 
tance in fields as diverse as optoelectronics and materials 
research, in biological applications such as spectroscopy 
and microscopy, in high harmonic generation, and in pho- 
toionization or molecular dissociation, among others. 

It is known that when the pulses duration approach 
the duration of only several optical cycles, theories based 
on the slowly-varying envelope approximation (SVEA) 
and the rotating-wave approximation (RWA) are not ap- 
propriate, since phenomena such as electric field time- 
derivatives leading to carrier- wave reshaping or the 
generation of high spectral components |9(, can not be 
described by such theories. In these situations, accurate 
numerical modeling based on the finite-difference time- 
domain (FDTD) method [3 is being increasingly used 
to investigate the full set of optical Maxwell-Bloch equa- 
tions HHHim. 

In this Letter, we investigate the possibility of phase 
control of transient four-wave mixing interactions in two- 
level systems driven by two mutually coherent intense ul- 
trashort pulses of central angular frequencies lu and 3a; . 
We employ a standard predictor-corrector FDTD numer- 
ical approach which incorporates all propagation effects 
-such as nonlinearity, dispersion, absorption, dissipation, 
saturation, and other resonance effects- || , to study the 
sensitivity on the relative phase (f> of the nonlinear cou- 
pling of the lu-3lu pulses to the field at frequency 5ui, 
which results from the interaction of the waves thr oug h 
the third order susceptibility of the medium (x^) |l2| . 

By describing the evolution of the field spectrum as 
the two-color pulses propagate through the medium, our 
simulations demonstrate that the generation of the Sub- 
component is enhanced when the two pulses are in phase 
((f> = 0), while it is nearly suppressed when they are out 
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FIG. 1: Schematic energy level diagrams corresponding to 
the first case of our study. The four-wave mixing can occur 
between the injected two-color fields at frequencies lu and 3lu, 
and a third generated field at frequency 5lu, by parametric 
coupling (a) or by two-photon processes (b). The wavelengths 
associated to the frequencies in this scheme are: Ar^) — 2400 



nm; A 



(3") 



800 nm; A 
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480 nm; A 



area of both pulses is A = 207r. 
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of phase ((f) — ir). In what follows, we analyze this effect 
in two different configurations: 

- First, we consider a nonresonant case with large pulse 
intensities (see Fig. |JJ. The two injected pulses have 
a duration of 10 fs and they initially overlap in time. 
The pulse area for both pulses is A = 207T |l3j, and the 



wavelengths are A 



O) 



2400 nm and A 



(3w) 



800 nm. 



The fields are largely detuned from the atomic transition 
resonance (Xt Ul2 -\ = 2-KcjuiYi = 400 nm), although the 
two-photon coupling of the field with = 800 nm 

is considered in resonance (see Fig. [I](b)). Therefore, 
by four-wave mixing, the fields can interact with a gen- 
erated third wave with A( 5w j = 480 nm, by either the 
parametric coupling represented in Fig. ^ ( a ) or through 
a two-photon process [Fig. 2] (b)] . We note that the en- 
ergy level diagrams of this first study can be considered 
as an example of the vibrational nonresonant contribu- 
tions at zero pump-probe time delay in CARS processes 

- In a second study, we consider two 10 fs pulses with 
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moderate pulse areas (A 



A 



1743 nm and A 



(3w) 



47r) and wavelengths as 
581 nm. In this case, the 
field with A( 3tJ ) = 581 nm propagates close to resonance 
with the atomic transition (A( Wl2 ) = 580 nm). By a non- 
linear third order parametric coupling (see Fig. |SJ) the 
fields can produce a wave with A( 5tlJ ) = 348.6 nm. The 
parameter values in this second study have been chosen 
to be close to those for transitions that are relevant in 
biological applications, with similar schemes being cur- 
rently investigated with femtosecond CARS techniques 

The Maxwell-Bloch equations can be written as |8| 
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where H y (z,t) and E x (z,t) represent the magnetic and 
electric fields propagating along the z direction, respec- 
tively, /io and eo are the magnetic permeability and elec- 
tric permittivity of free space, respectively, N at repre- 
sents the density of polarizable atoms, T is the dipole 
coupling coefficient, T\ is the excited-state lifetime, T2 is 
the dephasing time, io\2 is the transition resonance an- 
gular frequency of the two level medium, and p\ and p 2 
are the real and imaginary components of the polariza- 
tion -which determine the index of refraction and gain 
coefficients, respectively- @. The population difference 
is P3, and pzo represents its initial value. 

As mentioned above, an hyperbolic secant two-color 
pulse, which can be expressed as 

E x (z = 0,t)= E u (t) + E 3u) (t) = E sech({t - t )/t p )x 
[cos(uj{t - t j) + cos{3uj(t - t ) + 4>)\ , (2) 

is externally injected to the system. The peak input 
electric field amplitude Eo is chosen the same for both 
pulses and, in terms of the pulse area j4, it is given by 
E = Ah/{TrTt p ) [3. The duration of the pulse is given 
by t p = Tp/1.763, with r p being the full width at half 
maximum (FWHM) of the pulse intensity envelope, to 
gives the offset position of the pulse center at t = 0, and it 
is the reference value for the absolute phase of the pulses. 
The angular frequencies of the pulses are to and 3lu, and 
<j) is the relative phase. The parameters have been chosen 
as follows: T x = T 2 = 1 ns, r p = 10 fs, T = 2.65eA, and 
N at ~ 2 x 10 18 cm~ 3 . It is important to note that our 
simulations give information on the coherent (transient) 
behavior of the system, since the dephasing time scales 
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FIG. 2: Superposition of the electric field of the pulses. Eu(t) 
is the electric field of the pulse with frequency u represented 
as a function of time, and E^{t) represents the pulse with 
frequency 3w (see also Eq. (J5|. As indicated, the upper figures 
correspond to the addition of the pulses in phase ((/> — 0), 
and the lower ones correspond to the addition out of phase 

(4> = w). 



are chosen much larger than the duration of the pulses 
Ti,T 2 » t p . 

Figure [21 is a representation of the electric field of the 
injected pulses. In the two upper figures we can observe 
the addition of two pulses with frequencies u> and 3lu 
in the case that the relative phase 4> between them is 
4> = 0. We note that the constructive interference in the 
center of the pulses results in a maximum of the resulting 
pulse (Fig. [3 upper- right), while in the case that the 
pulses are added out of phase (</> = tt) , there is destructive 
interference in the center of the pulses and therefore the 
resulting pulse (Fig. [21 upper-right) presents a zero in 
the center. We will show below that the different shape, 
which only depends on the relative phase between the two 
lu — 3cj pulses, is crucial in the nonlinear interaction with 
the two-level medium in the case of ultrashort intense 
pulses. 

We next show the results of the simulations corre- 
sponding to the first scheme (Fig. ^) that we have ana- 
lyzed. In Fig. [21 we represent the spectrum of the field 
for different positions of the propagation. The two high- 
est peaks on the left of each spectra, in Fig. correspond 
to the frequencies of the pulses initially injected to the 
medium (u)/2tt and 3w/27r). When the two pulses are 
initially in phase (</> = 0), a third frequency (anti-Stokes) 
component is generated at 5w = 3w + (3w — w), as shown 
in Fig. Ela)- It is clear in that case that the conversion of 
the third harmonic component (3<x>) to the 5w-component 
[which is marked with a dotted arrow in Fig. Ofa)] in- 
creases as the propagation length increases. Differently, 
however, when the two pulses are initially out of phase 
(4> = 7r), the conversion to the 5cj-component is com- 
pletely suppressed, as it can be observed in Fig. Ofb). 
This contrasting behavior, to the best of our knowledge, 
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has not been observed before. We attribute this phe- 
nomenon to the dependence on the shape of the carriers 
(see Fig. [5J of the quantum interferences of the medium, 
which in the case of intense ultrashort pulses are governed 
by the electric field itself rather than by its envelope . 
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FIG. 4: Evolution of the field spectra as the pulses propagate 
through the two-level medium, (a) z — 50 /im; (b) z = 200 
fj,m; (c) z — 500 /im. Note that the verical axis is shown in 
logarithmic scale. 



FIG. 3: Evolution of the field spectra as the pulses propa- 
gate through the two-level medium, as indicated, (a) = 0; 
(b) <f> — tv. The dotted arrow in (a) indicates the peak at 
frequency 5u>. The spectra have been shifted right and up for 
different z. 



For longer propagation distances (z w 500 /im), the 
coupling between the fields will not longer be effective, 
because of changes in the shape of the pulses, as a conse- 
quence of the propagation effects, and also because of the 
progressive change in the spectral energy distribution of 
the field. In Fig. we have plotted the spectra of the 
propagating field, in logarithmic scale, for propagation 
distances up to z = 500 /im. We can see how the usual 
higher spectral components are increasingly generated in 
the field as the propagation distance increases. These 
higher spectral components are due, on the one hand, 
to the interference between the pulses and the associ- 
ated nonlinear processes 0]> and on the other hand, to 
carrier-wave Rabi flopping phenomena which are present 
due to the high energies contained in the injected pulses 
0. As can be observed in Fig. EJc), in the case of the 
present simulations, the spectrum spreads over the entire 
UV region for z = 500 /im. Note that some sensitivity on 
the initial relative phase of the spectra, in the higher 
frequency region, can also be observed in Fig. 0] Fur- 
ther details on these long propagation effects, however, 



are left to be given elsewhere. 




FIG. 5: Energy level scheme representing our second study. 
As in the previous scheme Q , the coupling can occur between 
the fields injected at frequencies u) and 3u, and a third field 
generated at 5w. The wavelengths associated to the frequen- 
cies in this scheme are: A(^) = 1743 ran; -^(3^) = 581 nm; 
-^(5u) = 348.6 nm; A(„ 12 ) = 580 nm. The area of the pulses 
is A = 4-7T. Note that, in this case, the field at 3a; is close to 
resonance with the transition. 



The second energy level configuration that we have an- 
alyzed is represented in Fig. [S] In this case, two overlap- 
ping pulses with a duration of t p = 10 fs are injected to 
the two-level medium. The area of the pulses is in that 
case much smaller (A = An) than in the previous scheme, 
and the field at frequency 3a; is now nearly in resonance 
with the transition of the medium. In Fig. the field 
spectra for two different positions are shown (z = 50 /im 
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[Fig. Eta)] and z = 300 [Fig. Et b )D- We clearly 
observe that the transfer of energy to the 5w-component 
(which is marked with a dotted arrow in Fig. EJ) is only 
efficient in the case of constructive interference in the cen- 
tre of the pulses (<fi = 0), while it is almost suppressed 
for (<fi = 7r). On the other hand, the rest of the main 
peaks in the spectra, which result from the interference 
of the fields, remain basically insensitive to the variations 
of <j), and we therefore conclude that the influence of the 
initial relative phase is mostly observed in the four-wave 
coupling to the 5w-component. This second study hence 
provides another demonstration of the possibilities of the 
phase control phenomena that we report. 
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FIG. 6: Field spectra at z — 50 /im (a) and z = 300 /im 
(b). The dotted arrows indicate the peaks at frequency 5a; 
for <j> = 0. 

We finally note that the results of the second scheme 
that we have investigated can be compared with previous 
work , where in a similar configuration the production 
of higher spectral components even for small pulse areas 
was demonstrated, and it was shown that the oscillatory 
structures around the resonant frequency and the prop- 
agation features of the laser pulses depend sensitively on 
the relative phase <f> of the two pulses. We have here 
addressed the influence of the relative phase on the non- 
linear coupling to the optical field at frequency 5u> , which 
was not considered in [7j, and in that sense the second 
study in our work can be considered as complementary 
to that investigation. 

In conclusion, we have analyzed the coherent propaga- 
tion of two-color ujSlj femtosecond laser pulses overlap- 
ping in time and propagating in a two-level system to an 
extend of some hundreds of microns. Our study predicts 



a critical dependence on the relative initial phase for the 
transfer of energy by four-wave mixing to the field at 
frequency 5uj. This effect is clearly observed in our sim- 
ulations for propagation distances as long as z « 300 
/im. We have observed this phenomenon in the case of 
highly intense femtosecond pulses strongly detuned from 
the transition of the medium, and also for pulses with 
moderate intensity close to resonance with the material 
transition. We have hence demonstrated that the manip- 
ulation of the initial absolute phase of the different fre- 
quency components of sinchronously propagating uj-3uj 
femtosecond pulses can serve as a means to control the 
nonlinear coupling to the optical field at frequency bco in 
four-wave mixing interactions. This coherent control can 
be useful, in particular, to limit the background nonres- 
onant contributions in coherent anti-Stokes Raman pro- 
cesses. 
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